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Abstract 

The cross section of neutron-deuteron radiative capture nd ^Hj is calculated 
at energies relevant to Big-Bang nucleosynthesis { 20 < E < 200 KeV ) with pi- 
onless Effective Field Theory. At these energies, magnetic transition Ml gives the 
dominant contribution. The Ml amplitude is calculated up to next-to-next-to leading 
order(N^LO) with insertion of three-body force. Results are in good agreement within 
few percent theoretical uncertainty in comparison with available calculated data below 
E=200 KeV. 
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1 Introduction 



Very low-energy radiative capture and weak capture reactions involving few-nucleon systems 
have considerable astrophysical relevance for studies of stellar structure evolution and big- 
bang nucleosynthesis. At these energies pionless Effective field theory(EFT) is an important 
tool for computing physical quantities. 

Much of the strength of EFT lies in the fact that it can be applied without off shell ambi- 
guities to systems with more nucleons. In past years, nuclear EFT has been applied to two-, 
three-, and four-nucleon systems[l-9] and recently developed pionless EFT is particularly 
suited to high-order and precision calculation. An example of a precise calculation is the 
reaction np — >■ 'yd, which is relevant to big-bang nucleosynthesis(BBN), and the cross section 
for this process was computed to 1% error for center of mass energies E <1 MeV jlOj . 

On the other hand. The three-body EFT calculations have been so far confined to 
nucleon-deuteron system. For example nucleon-deuteron scattering in all channels expect 
the 5*1/2- wave can be calculated to high-orders using only two-nucleon input and the triton 
binding energy is found to be B^^^^'^^ = 8.35MeV in next to leading order close to the ex- 
perimental Bs^^'^P^ = 8.5MeV [Hj. Pionless EFT has been recently applied to the four-body 
system with contact interactions and large scattering length at leading order by Platter [T^ . 
The binding energies of the ^He tetramer and the alpha-particle have been calculated and 
it is in a good agreement with experimental value. 

The radiative capture of neutrons and protons by deuterons and the inverse reactions, the 
photodisintegration of and ^He, have been investigated experimentally and theoretically 
over the last decades with some interest. In an experiment performed in recent years at 
TUNL ^21 total cross section and vector and tensor analyzing power of the pd — 

^He^ process were measured at the center of mass energies below 55 KeV. 

The theory of the nd — > ^H'y capture reaction has a long history. The nd doublet 
or quartet state Ml transition calculated in Impulse Approximation, and explanation of 
smallness of its cross section when compared to the np —>■ dj reaction, ax = 334.5 ±0.5 mb, 
were pointed out by Schiff ^3]. Later, Phillips fH] emphasized the importance of initial state 
interactions and two-body currents to capture reaction in the three-body model calculation, 
by considering a central, separable interaction. In recent years, a series of calculations of 
increasing sophistication with regard to the description of both the initial and final state 
wave functions and two-body current model were carried out fTj . These efforts culminated 
in the 1990 Friar et al. jTH] calculation of the nd ^H-y total cross section, quartet capture 
fraction, and photon polarization, based on converged bound and continuum state Faddeev 
wave functions, corresponding to a variety of realistic Hamiltonian models with two- and 
three-nucleon interactions, and a nuclear electromagnetic current operator, including the 
long-range two-body components associated with pion exchange and the virtual excitation 
of intermediate A resonances. 

For very low energy the p-d and n-d radiative capture, a magnetic dipole(Ml)transition 
is a dominant contribution,which was studied in plane wave (Born) approximation by Friar et 
al. ^H] • In these investigations the authors employed their configuration-space Faddeev cal- 
culations of the helium wave function, with inclusion of three-body forces and pion exchange 
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currents. Various trends, e.g., the correlation between cross sections and helium binding en- 
ergies, and their potential dependence were pointed out. More recently a rather detailed 
investigation of such processes has been performed by Viviani et al. UHl- Their calculations 
employed the quite accurate three-nucleon bound- and continuum states obtained in the 
variational pair-correlated hyperspherical method, developed, tested and applied over years 
by this group. 

We calculate very low energy cross section of radiative capture of neutrons by deuterons 
nd — i> 7^ if at energies 20 < < 200 KeV, relevant to Big-Bang nucleosynthesis, with 
pionless EFT. At these energies, magnetic transition Ml gives the dominant contribution. 
The Ml amplitude is calculated up to next-to-next-to leading order(N^LO) with insertion of 
three-body force. Results show good agreement in comparison to ENDF [20] below E=200 
KeV. 

The organization of the paper is as follows: We first describe the relevant Lagrangian 
and scattering in doublet S-wave channel in Section |21 This section essentially introduced to 
define various parameters that enter in the expression for the cross section. The calculation 
of the total cross section is presented in sections El We tabulate the calculated cross sections 
for some energies relevant for BBN, discuss the theoretical errors, and compare our results 
with the corresponding values from the on-line ENDF/B-VI database [201 Section HI 
Summary and conclusions follow in Sectional 

2 ^/S'i/2 neutron-deuteron scattering(triton channel) 

The ^Si channel to which ^He and belong is qualitatively different from the other three- 
nucleon channels. Consequently, describes the preferred mode for nd — ^H'-y and 
pd — >■ ^He'y processes. This difference can be traced back to the effect of the exclusion 
principle and the angular momentum repulsion barrier. In all the other channels, it is either 
the Pauli principle or an angular momentum barrier(or both)which forbids the three-particle 
to occupy the same point in space. As a consequence, the kernel describing the interaction 
among the three-nucleon, unlike in the bosonic case, is repulsive in these channels. The zero 
mode of the bosonic case, describes a bound state since it is a solution of the homogeneous 
version of the Faddeev equation. As such, it is not expected to appear in the case of repulsive 
kernels and, in fact, it does not. For the ^He or channel however, the kernel is attractive 
and, as we will see below, closely related to the one in the bosonic case [TT] . 

Let us first discuss the integral equation describing nucleon-deuteron scattering. We 
start with the three-nucleon lagrangian which is given by |llj : 

C = Art + ^) AT + df (-i9o -J^ + ^s)dU d^ (-i9o - ^ + A.) dl 

+ ('^° + ^ + M + ^ " (rff (AT^P^AT) + H.c.) - gt (dfiN^'P^N) + H.c.) 
- uj, {t\T^N)dj + H.c.) - ut {t\a'N)di + H.c.) + . . . , 

(2.1) 
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where is the nucleoli iso-doublet and the auxihary fields t, df and dl carry the quan- 
tum numbers of the '^He-^H spin and isospin doublet, ^5*0 di-nucleon and the deuteron, 
respectively. The projectors P* and are defined by: 

P^ = -^^2r2T^ , (2.2) 

V o 

where A = 1,2,3 and i = 1, 2, 3 are iso-triplet and vector indices and (a*) are isospin 
(spin) Pauli matrices. 

One can write the Faddeev integral equation in the kinematics defined by the two 
cluster-configurations exist in the three-nucleon system. We follow the notation suggested 
by GrieBhammer in [21]. The A^cij-cluster with total spin S* = | or = |, depending on 
whether the deuteron and nucleon spins are parallel or anti-parallel; and the A^rf^-cluster 
which has total spin S — ^, as d'^ is a scalar. The leading-order three-particle amplitude is 
0{Q~^) (before wave-function renormalisation) and includes all diagrams built out of the 
leading two-body interactions. The resultant Faddeev integral equation is represented in 
Fig. As the Lagrangian up to N^LO does not mix partial waves or fiip the spin of the 




Figure 1: The Faddeev equation for A^d-scattering up to N^LO. Thick solid line: propagator 
of the two intermediate auxiliary fields ds and dt, denoted by V, see ()2.5p : /C: propagator 
of the exchanged nucleon; 7i: three-body force and the cross denotes insertion of deuteron 
kinetic energy operator. 

auxiliary fields, angular momentum is conserved in the quartet and doublet channels. 

In the doublet channel, the Faddeev equation is two-dimensional in cluster-configuration 
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space as both Ndt- and A^c/g-configurations contribute |21j : 
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is built out of the two amphtudes which get mixed: td^tt for the Ndt 
td^ts for the Ndt A^c/^-process. Furthermore, 



Dt{po,p) 








(2.4) 

A^rft -pro cess, and 
(2.5) 



is the propagator of the two intermediate auxihary fields. 

One included the specific three-body force term in the doublet-channel Faddeev equation 
(|2.3p at a given order n of expansion as suggested by Bedaque et al. jH 13 HTJ if and 
only if that term is needed to cancel cut-off dependences in the observables which are 
stronger than cut-off dependence from the suppressed terms of order {QR)^~^^ where Q is 
a typical external momentum and R is the short distance scale in position space beyond 
which the EFT breaks down. This argument suggests that from the point of view of the 
EFT, cutoff dependences cancel order by order in the expansion. One can implement this 
idea by expanding the kernel of the integral equation perturbatively, and then iterate it 
by inserting it into integral equation. This partial re-summation arbitrarily includes higher 
order diagrams and consequently dose not improve the precision of calculation. The only 
necessary re-summation is the one present at LO. The partial re-summation of range effects 
is made only for convenience |3]. 

The integral equation is solved numerically by imposing a cut-off A. In that case, a 
unique solution exists in the ^5'i/2-channel for each A and H = 0, but no unique limit as 
A — > cxo. As long-distance phenomena must however be insensitive to details of the short- 
distance physics (and in particular of the regulator chosen), Bedaque et al. jH CH 121] 
showed that the system must be stabilized by a three-body force 



n=0 



(ME+jI 
V A2 



2ifo(A) , 2H2iA) 



A2 



A4 



(ME + ^l) + 



(2.6) 



which absorbs all dependence on the cut-off as A ^ cxd. It is analytical in E and can be 
obtained from a three-body Lagrangian, employing a three-nucleon auxiliary field analogous 
to the treatment of the two-nucleon channels [TT] . 
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H2 is dimension-less but depends on the cut-off A in a non-trivial way, as a renormali- 
sation group analysis reveals: Instead of approaching a fixed-point as A oo, it shows an 
oscillatory behavior known as "limit cycle" 0. 

As one needs a three-body force at LO, Hq ~ Q~^, all three-body forces obtained by 
expanding Ti. in powers of E are also enhanced, with the interactions proportional to H2 
entering at N^LO fJE]- The power-counting for the three-body forces is hence 

Ho{A)^Q~' , H2{A)^Q-' , i/2n(A)~g-^ (2.7) 

The scattering phase-shift of the S-wave in the quartet and doublet channel is related 
to the renormalised on-shell amplitudes by 

Sir 1 Stt 1 

Tg = Zttq = — ^ , T,,,y = — -^—^——^ , (2.8) 

where x,y = s,t label the matrix entries in cluster-configuration space, and 

fd = ZU with Z := 1^ , (2.9) 



,0 ^%Z, 

is the renormalised doublet-amplitude and its wave-function renormalisation. In the doublet 
channel, the only observable process is nucleon-deuteron scattering, Ndt Ndt, i.e. x = 
y = t. 

Nucleon-deuteron scattering is to N^LO thus completely determined by four simple ob- 
servables of A^A^-scattering: the deuteron binding energy,residue (or effective range), the 
scattering length and effective range of the ^5'o-channel. Only the ^S'i/2-channel has further 
unknowns, namely the strength of the three-body interaction Hq at LO and NLO, and in 
addition of H2 at N^LO . They are determined by its measured scattering length ad |22] and 
the triton binding energy B^, respectively: 

ttd = (0.65 ± 0.04) f m , = 8.48 MeV . (2.10) 

3 Neutron-deuteron radiative capture process 

In the KeV energy region, the processes nd — s> ^H'j and pd — > ^He'y play an important role 
in nuclear astrophysics and in nuclear physics. In the standard big-bang nucleosynthesis 
theory the corresponding reaction rates are necessary to estimate the ^He-yield as well as 
the abundances of other light elements. The spin structure of the matrix elements of neu- 
tron radiative capture by deuteron is complicated but in very low energy for this reaction 
we can introduced three multipole transition that can be allowed by p-parity and angular 

-J _|_ Q -|- 

momentum conservation i.e. = ^ — Ml and = | — Ml, £'2. The parameteriza- 
tion of the corresponding contribution to the matrix element to be build by the following 
contributions: 

i{t^N){D ■ e* X k) , 
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{t^aaN){D X [e* x ^]), , (3.1) 

t\a-e* D-k + a-k D-e')N , 

where N, t, e, D and k are the 2-component spinors of initial nucleon filed, final '^He 
(or ^i7)field, the 3-vector polarization of the produced photon, the 3-vector polarization of 
deuteron and the unit vector along the 3-momentum of the photons, respectively. The two 
structures in Eq. p.ip correspond to the Ml transition. At thermal energies the reaction pro- 
ceeds through S-wave capture predominantly via magnetic dipole transition, M/"^*^, where 
L=0, S'=l/2,3/2 and i=l. To obtain the spin structure, which corresponds to a definite 
value of J for the entrance channel, it is necessary to build special linear combinations of 

- ^ - 1+ 3 + 

products DN and a x DN, with ~ 2 or = - : 

01/2 = {iD + ax D)N and {2iD -ax D)N . 

1+ 3 + 

For both possible magnetic dipole transitions with = - (amplitude gi) and = - 

(amplitude gs) we can write: 

gi : t\iD ■e*xk + axD-e*x k)N, 



gs: t\iD-e* xk + ax D-e* xk)N . (3.2) 

The electric transition E^^^ for energies of less than 60 KeV dose not contribute to 
the total cross section. Therefore transition will not be considered in energies 

relevant to BBN calculation. The Ml amplitude receives contributions from the magnetic 
moments of the nucleon and dibaryon operators coupling to the magnetic field, which are 
described by the lagrange density involving fields: 

= ^^^(^0 + ^^^'^-^ + ^7:r^4^=TWT^*'^^«35, + H.C . (3.3) 

where the unknown coefficient Li, which contributes at order Q must either be predicted 
from QCD or determined experimentally in order to have model-independent predictive 
power. 

The radiative capture cross section nd ^H'-f at very low energy is given by 



where 



PfiN 
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Figure 2: The Faddeev equation for A^(i-radiative capture up to LO. Boxes indicate insertion 
of A^ci-scattering amplitude up to LO from Fig. CI only up to first line of perturbative ex- 
pansion of the Faddeev equation). Cross shows wave function renormalization in each step 
when triton is made and three-body interactions are shown with strength Ho{A). Wavy line 
shows photon and small circles show magnetic photon interaction. The photon is minimally 
coupled. Remaining notation as in Fig. ^ 
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Figure 3: The NLO contribution can then be obtained by perturbing around the LO 
solution with the one deuteron kinetic energy operator insertion and last lines show NLO 
order of photon interaction with the lagrangian Eq. ()3.3|) with the Li vertices. Boxes indicate 
insertion of iV(i-scattering amplitude for NLO from Fig. ^only second line of perturbative 
expansion of the Faddeev equation). Remaining notation as in Fig. |21 

with X stands for either the electric or magnetic transition and fijsf is in nuclear magneton 
and p is momentum of the incident neutron in the center of mass. There is an infinite number 
of diagrams contributing at leading order for Ml amplitude of the radiative capture cross 
section nd — *• ^Hj, as shown in Fig. |21 We follow the same procedure as First, expand 
the kernel of the integral equation perturbatively in the same way of Fig. ^ and then iterate 
it by inserting it into the integral equation. We use the deuteron propagator in the integral 
equation, which is then solved numerically. The re-summation is necessary since according 
to the power counting all these diagrams contribute equally to the final amplitude. The 
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Figure 4: The N^LO contribution can then be obtained by perturbing around the LO 
solution with the twice deuteron kinetic energy operator insertion. Boxes indicate insertion 
of A^(i-scattering amphtude to N^LO from Fig. ^ Remaining notation as in Fig. |21 and 
Fig. El 




( ^' ) ( b ) 

Figure 5: The two different contribution of N^LO with twice deuteron kinetic energy oper- 
ator insertion. 

NLO and N^LO contributions can then be obtained by perturbing around the LO solution 
with the one and twice deuteron kinetic energy operator insertion, respectively, as shown 
diagrammatically in Fig. IHl and Fig. |31 For N^LO contribution due to very cumbersome 
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numerical calculation, one must substitute diagram of Fig. Efa) by Fig. Efb) in diagrams 
of Fig. I3J particularly when one needs to compute the full off-shell LO amplitude before 
inserting the deuteron kinetic energy terms. 

We now turn to the amplitude to be used in the Ml calculation. The amplitude td 
after properly iterated will be folded to electromagnetic interaction JCem order by order and 
integrated on momentum in each order. Finally the wave function renormalization in each 
order will be done when triton is made, we introduce for JCem '■ 



- rri's.) --pL_ + l,(.s,)|q|2 _ ,|q| 



(3.6) 



We have not considered the two following possible diagrams in our calculation. Photo- 
interaction directly with exchanged nucleon is shown in Fig. |Bfa). For this interaction, we 
have p. A and in very low energy relevant to BBN, p < 200 KeV ,this particular contribution 
(Fig. Efa)) which could be appeared in Figures (2,3,4) has been neglected, because we can 
estimate it's size from a naive power counting. In the integration of photon coupling to 
exchanged-nucleon vertices, we use q ^ Q and q ^ p, where we have the Q/M^y'^ the wave 
function renormalization, M/Q'^ nucleon propagator and l/My^Q dibaryon propagator. 
Consequently, it's size of contribution has been estimated to be < 1%. The other direct 
interaction is with three-body vertices H, up to order of our calculation only Hq, H2 are 
considered. The direct interaction in order of Hq, for zero energy range, will be dealt with 
elsewhere. This diagram is shown in Fig. lEfb). Contribution of this diagram for the energy 
range of our calculation is ignored, because in very low energy for p ~ g ~ Q error for 
neglecting the contribution of this term is less than 1%. Contribution of diagrams with H2 
vertices generally appears at N^LO. At N^LO one can insert photon also to H2 vertices. 
From ()2.6|1 . this contribution is (p^/A^)if2 and is of the order of {k/AY or (Q/A)^ by power 
counting expansion. In numerical calculation we define H2 such that it dose not contribute 
at zero momentum, naturally for the energy range near zero momentum, insertion of photon 
to H2 vertices for momentum p ~ 200 KeV, could be neglected and the error is estimated 
to be smaller than the estimation of error for insertion of photon to Hq vertices. 

The solution depends on A for arbitrary Hq, H2 which in general are cutoff dependent, 
too. Since the low-energy amplitude can not depend on the cutoff, any explicit dependence 
on the cutoff has to be canceled by the cutoff dependence implicit in Hq, H2, etc. We 
determine the cutoff dependence of the three-body forces by imposing this condition order 
by order in 1/ A in the same way of . In this process, we also determine which three-body 
force appears at every order of the expansion. 

In order to have cutoff independent amplitude one must determine the three-body forces 
by demanding order by order in Q/A 



A 

,2 



n{E,A) + ^ j dq q^V{E - ^,q)[IC{q,p) +n{E,A)]t^{q) = const. , (3.7) 
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(■•0 (b) 

Figure 6: Diagram (a) shows photon interaction with exchanged nucleon and diagram (b) 
shows interaction of photon with Hq vertex. 



This (possibly energy dependent) constant can without loss of generality be set to zero by 
absorbing it into a re-defined three body force 



const. 



1 + ^ j dq q^V{E 



(3.8) 
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Figure 7: Semi-logarithm curve of the cutoff variation of phase shifts is shown between 
A = 150 MeV and A = 500 MeV as function of the center-of-mass momentum . The short 
dashed, long dashed and solid line correspond to LO, NLO and N^LO, respectively. 



Imposing ()3.7|) with const. = order by order in the low energy expansion, we can 
determine which three-body forces are required at any given order, and how they depend 
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on the cutoff. As we are interested only in the UV behavior, i.e. in the asymptotics, and 
no IR divergences occur, we can safely neglect the IR limit of the integral. All calculations 
were performed with the cutoff at 400 MeV. At N^LO , where we saw that H2 is required 
the cutoff is varied between 150 and 500 MeV and we ignore other sources of error at low 
cutoffs, this is a reasonable estimate of the errors. For energy range of our calculation, 
below the break up reaction, we show cutoff variation of phase shifts between A = 150 
MeV and A = 500 MeV as function of the center-of-mass momentum in Fig. [7| We can see 
very smooth slope and does however change significantly from order to order because the 
dominant correction is (7/A)" and for these energy range at every order these variations are 
nearly independent of variation of momentum. We confirm that in our calculation also in 
very low energy range the cutoff variation decreases steadily as we increase the order of the 
calculation and it is of the order of (/c/A)", (7/A)", where n is the order of the calculation 
and A = 150 MeV is the smallest cutoff used. It is worth mentioning that the errors due to 
increasing momentum, as one would expect due to (A;/ A)" also appear in our calculation in 
very low energy but these errors decrease when the order of calculation is increased up to 
N^LO. 

4 Numerical results for neutron-deuteron radiative cap- 
ture 

We numerically solved the Faddeev integral equation up to N^LO. We used fic = 197.327 MeV fm, 
a nucleon mass of M = 938.918 MeV, for the A^A^ triplet channel a deuteron binding energy 
(momentum) of B = 2.225 MeV (7^ = 45.7066 MeV), a residue of = 1.690(3), effective 
range rot = 2.73 fm, for the A^A^ singlet channel an ^Sq scattering length of at = —23.714 fm 
and fiN = 5.050 x 10"^^JT~^. We determine the two-nucleon parameters from the deuteron 
binding energy, triplet effective range (defined by an expansion around the deuteron pole, 
not at zero momentum), the singlet scattering length, effective range (defined by expanding 
at zero momentum), and two body capture process(obtained with comparison between ex- 
perimental data and theoretical results for nd — d'-f process at zero energy jTHI). We fix the 
three-body parameters as follows: because we defined H2 such that it does not contribute 
at zero momentum scattering, one can first determine Hq from the ^Si scattering length 

as = (0.65 ±0.04) fm 122] • 

At LO and NLO, this is the only three-body force entering, but at N^LO, where we saw 
that H2 is required, it is determined by the triton binding energy B-^ = 8.48 MeV. We 
solve integral equation by expansion in order of Q and properly iterating the kernel. Then, 
the resulted will be folded to electromagnetic interaction order by order and properly 
integrated on the involving momentum. 

For our calculation source of error due to low cutoffs and low momentums for very low 
energy calculation is neglected. The cutoff variation decreases steadily as we increase the 
order of the calculation and is of the order of (fc/A)", (7/A)", where n is the order of the 
calculation. We used A = 150 MeV, for the smallest cutoff . 

The cross section calculation for neutron radiative capture by deuteron as function of 



12 



ID 



s 




LO 








NLO 








N^LO 




6 
4 








2 

























05 075 0.1 0.125 0.15 175 3 

E[MeV] 



Figure 8: The cross section for neutron radiative capture by deuteron as function of the 
center-of-mass kinetic energy E in Mev. The short dashed, long dashed and sohd hne 
correspond to the contribution from Ml capture up to LO, NLO and N^LO, respectively. 



Table 1: Neutron radiative capture by deuteron in micro barn and estimated error in per- 
centage in comparison with the last column at every order up to N^LO. Last column shows 
ENDF results for cross section |2Uj . 



Energy 
(KeV) 


LO 


error (%) 
LO 


a{fib) 
NLO 


error(%) 
NLO 


a{fib) 
N^LO 


error (%) 
N^LO 


ENDF{fxb) 


40 


1.64(6) 


30 


1.31(5) 


4 


1.25(0) 


1.5 


1.27(0) 


50 


1.72(8) 


24 


1.45(8) 


4 


1.38(3) 


0.7 


1.39(0) 


60 


1.90(5) 


27 


1.58(5) 


5 


1.50(1) 


0.1 


1.50(0) 


70 


2.02(6) 


25 


1.72(0) 


7 


1.61(6) 


0.1 


1.61(0) 


80 


2.11(9) 


22 


1.82(2) 


6 


1.73 (3) 


0.5 


1.72(0) 


100 


2.42(6) 


25 


2.04(4) 


5 


1.93(3) 


0.5 


1.94(0) 


140 


2.88(9) 


30 


2.41(8) 


9 


2.30(3) 


0.1 


2.22(9) 



the center-of-mass energy at LO, NLO and N^LO is shown in Fig.|Hl 

Table 1 shows numerical results for the EFT nd — ^ ^H'y cross section for various nucleon 
center of mass energies E up to N^LO and errors estimate at every order in comparison with 
the last column. The corresponding values for the cross section from the online evaluated 
nuclear data file ENDF/B-VI [20] are shown in the last column. The EFT results for this 
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Figure 9: Estimation of error in percentage is shown in comparison with ENDF [20] versus 
nucleon center-of-mass kinetic energy E in KeV. The short dashed, long dashed and sohd 
hne correspond to the error up to LO, NLO and N^LO, respectively. 

cross section are presented up to only two significant digits(the third digit is shown for better 
comparison with ENDF data). In Fig. IHl we show the error due to available evaluated data 
ENDF [20] in percentage versus nucleon center-of-mass kinetic energy E in KeV. 

5 Conclusion 

We have calculated the cross section of radiative capture process nd — > ^H'y. We applied 
pionless EFT to find numerical results for the Ml contributions for this capture process 
for incident neutron energies relevant for BBN, 0.02 < E < 0.2 MeV. At these energy our 
calculation is dominated by S-wave state and magnetic transition Ml contribution. 

The error estimate in the cross section in comparison with evaluated nuclear data file 
ENDF 120] is shown in Fig. Eland Table 1. Errors estimate are 20-30 percent at leading 
order, below 10 percent up to NLO and by insertion of three-body force at N^LO, this 
error is reduced to below 1% percent. Specially, our calculation shows minimum error for 
energy 60-70 KeV up to at N^LO. Comparison of the LO with the NLO and N^LO results 
demonstrate convergence of the effective field theory. Finally, three-body forces will enter 
at higher orders of the EFT approach and reduce the theoretical uncertainty. 
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